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The Hubble diagram extended to z » 1: the 7— ray properties of 
GRBs confirm the ACDM model 
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ABSTRACT 

Tight constraints on cosmological parameters can be obtained with standard candles spanning 
a range of redshifts as large as possible. We propose to treat SNIa and long Gamma-Ray 
Bursts (GRBs) as a single class of candles. Taking advantage of the recent release of the Su- 
pernova Legacy Survey and the recent finding of a tight correlation among the energetics and 
other prompt j~ray emission properties of GRBs, we are able to standardize the luminosi- 
ties/energetics of both classes of objects. In this way we can jointly use GRB and SNIa as 
cosmological probes to constrain f2 m and J7a and the parameters of the Dark Energy equation 
of state through the same Bayesian method that we have, so far, applied to GRBs alone. De- 
spite the large disparity in number (115 SN la versus 19 GRBs) we show that the constraints 
on D, m and Ha are greatly improved by the inclusion of GRBs. More importantly, the result 
of the combined sample is in excellent agreement with the ACDM concordance cosmological 
model and does not require an evolving equation of state for the Dark Energy. 
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1 INTRODUCTION 

Cosmology passed from being mostly a theoretical science to be 
one of the most accurate physical sciences in the phenomenologi- 
cal sense. The recent high-precision measurements of cosmologi- 
cal parameters together with the spectacular advances in the under- 
standing of cosmic structure formation, produced a coherent picture 
of the evolution of the Universe but, on the other hand, prompted 
new fundamental questions. One of them is related to the expansion 
history of the Universe and the possibility that a repulsive medium 
with negative pressure (Dark Energy, hereafter DE) dominates its 
content. The Hubble diagram (distance-redshift relation, HD here- 
after) provides a key probe of the Universe expansion history. Be- 
cause of their almost homogeneous intrinsic energetics, Type la Su- 
pernovae (SNIa) have been used as the main distance indicators for 
constructing the HD. However, their current observability is limited 
to redshifts z < 1.7, and the high-z SNIa could suffer intergalac- 
tic du st extinction and evolutive effects fe.g.. lOstman & Mortselll 
l2005t and the references therein). Thus, the finding of alternative 
and complementary cosmological distance indicators is highly de- 
sirable. Gamma-ray Bursts (GRBs), after "standardizing" their en- 
ergetics through adequate relations, probed to be reliable distance 
indicators detectable up to very high z's and free from dust ex- 
tinction. The recent discovery of a tight correlation among prompt 
7— ray emission observables alone iFirmani et alj|20 06ai opens 
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new perspectives for the HD method, where G RBs can be com - 
bined with SNIa to extend the HD to z ~ 5 - 10 JSchaefeJ2003Ti . 

The observations of SNIa demonstrated th at the exp ansion 
of the Universe is accelerating iRiess et allll998t IPerlmutter et alJ 
ll99St) . The main explanation to this acceleration is the dominion of 
DE in the current cosmological dynamics, though departures from 
conventional physics, like modified gravity theories, are also con- 
sidered. DE is characterized mainly by its equation-of-state param- 
eter, w — Pdb/pdb, where pde and pdb are the pressure and en- 
ergy density of DE. For w < —1/3 the universe undergoes accel- 
erated expansion. The simplest interpretation of DE is the cosmo- 
logical constant A, in which case w = — 1 and pde = PA =const. 
However, more exotic models, with w 7^ — 1 and in general varying 
with 2, have been proposed (e.g., quintessence, k-essence, Chaply- 
gin gas, Braneworld models, etc.); even models that allow w < — 1 
(e.g., Phantom energy ) have been con sidered (see for recent re- 
views ISahnl2004l;|Padmanabharl20o'6l) . 

We need further observations to unveil the true nature 
of what we call DE. A strong effort is being done now 
for developing t he next-generation SNIa experiments (see e.g. 
iLinder & HutereJ2005ti aimed mainly to reduce random uncertain- 
ties. However, it is also crucial to reduce systematic uncertainties 
as well as to break model degeneracies (e.g., IWejlex&^Al brechj 
2002; Linder & Huterer 2003; Nesseris & Perivolaroooulos 2005; 
Ghise llini et al J I2005 : Firm ani et all200 6hi. The two latter papers 
illustrate how some degeneracies in the cosmological parameter 
space can be reduced if the standard candles used to construct 
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the HD span a wide z range. In this sense, long GRBs have been 
proposed as a class of objects abl e to exte nd the HD up to very 
high redshifts in a natural manner (Schaefer 2003; Ghirlanda et alJ 
l2004bh . 

Despite the large disp ersion of the long G RB energetics 
iFrail et al J l200ll: iBloom. Frail & Kulkarni 2003), the discovery 
that their energetics correlate with observable quantities like the 
peak energy (in vL v ) of the time-integrated prompt emis- 
sion spectrum an d the achromatic "bre a k-time" in the after- 
glow ligth curve jGhirlanda etafl l2004at Liang & Zhand l2006t 
iNavaet alJ 2006), has been used to "standardize" it. This al- 
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Firmani et al. (2006a) have found a new 



GRB correlation whose tightness, in the framework of the standard 
fireball scenario, is explained by its scalar nature. Since such new 
correlation has the same shape in the observer and in the comoving 
frame, the influence of the T relativistic factor on the observed scat- 
ter becomes negligible. The correlation is based on prompt 7— ray 
observables only (besides the redshift), by-passing the need of 
measuring afterglow quan tities as is the case of the "Ghirlanda" 
relation ( Ghirla nda et alJ |2Q04a). The quantities involved in the 
new correlation are the isotropic peak luminosity Li so , _E p k, and 
the "high signal" timescale To. 45, previously used to c haracterize 
the variability behavior of bursts. In Firmani et al. (2006b) we have 
found that by varing the cosmology, the data points present a mini- 
mum scatter around their best fit line in correspondence of the so- 
called A cold dark matter (ACDM) concordance cosmology: a flat 
geometry Friedmann-Robertson-Walker-Lemaitre model with the 
cosmological constant dominating today. This result shows that, 
indeed, the I/i S o--E p k-To.45 relation can be used to derive cosmo- 
logical constraints. 

Due to the lack of low— z GRBs for calibrating a given cor- 
relation independently from cosmology, the use of a statistical ap- 
proach to jointly calibrate the correlation and c onstra in the cosmo- 
logical parameters is required. Firmani et al] |2005, 2006b) have 
presented an iterative Bayesian method to deal with this, so-called, 
'circularity problem'. The same method can be used for SNIa. 
Note that SNIa are not perfect standard candles: their luminosities 
vary with the shape of the light-curve (the brighter-slower rela- 
tion) and with color (the brighter-bluer relation) (e.g.. lGuv et alJ 
2005, and the references therein). Due to several high— z system- 
atical effects, a better calibration of these relations is obtained if 
higher— z SNe are included. The latter makes these relations cos- 
mology dependent. Therefore, the best fit to these rel ations has to 
be ca rried out jointly within the same cosmological fit lAstier et alJ 
I2005L hereafter A05). This approach has been applied to the 'Su- 
pernova Legacy Survey' (SNLS) of 115 SNIa (A05). We apply 
here our Bayesian method to this sample to improve the constraints 
given by A05, who used a simple multi-parametric \ 2 minimiza- 
tion method. 

Observations of SNIa are accurate and the current sam- 
ples comprise more than one hundred objects, but they are de- 
tected only at relatively low z's, which introduces the degen- 
eracy problem mentioned above. The GRBs useful as distance 
indicators span a large redshift range (up to z — 4.5) but 
they are still scarce (19 bursts for the I/i S o--E P k-Tb.45 relation). 
Thus, a promising strategy to partially overcome the problems 
that each family of objects individually suffer of, is to com- 
bine both in the same Hubble diagram and use them jointly for 
constraining the cosmological parameters. This is the goal of 
this Letter. Our main result is that the concordance A cosmol- 
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ogy (minimal DE model) is fully consistent with the joint GRB 
and SNLS SNIa data spanning the redshift range from z = 
0.015 to 4.5. Previous results with the so-called SNIa 'gold-set' 
(z < 1.7) showed a marginal inconsistency with the concor- 
dance model (Riess et al. 2004; Alam, Sahni, & Starobinskv 200' 
Choudhurv & Padmanabhan 2004; Jassal, Baela & Padmanabha: 
2005; Nesseris & Perivolaroooulos 2006), suggesting the possibil- 
ity of alternative DE models. 

In §2 we describe the SNIa and GRB samples used here. Sec- 
tion 3 deals with the method while the results are presented in §4 
The conclusions of this work are given in §5. 



2 THE SAMPLES 
2.1 Type la supernovae 

We analyze here the SNIa sample presented in A05. This sample 
consists of 44 nearby (0.015 < z < 0.125) SNIa assembled from 
the literature, and of 73 distant SNIa (0.249 < z < 1.010) dis- 
covered and followed during the first year of SNLS 1 . The same 
light-curve fit method (Guv et al. 2005) was applied to all SNIa in 
the sample. For each SN, the reported quantities along with their 
errors are the fitted rest-frame B— band magnitude m* B and the 
values of the parameters s (light-curve stretch) and c (normalized 
B — V color at the maximum of the light-curve). The magnitude 
m* B refers to observed brightness, and therefore does not account 
for the brighter-slower and the brighter-bluer correlations. As men- 
tioned in the Introduction, A05 suggest to empirically calibrate 
these correlations using all objects (either at low or high— z). For 
the cosmological fits, two SNIa out of the 117 SNLS sample objects 
were excluded because they are outliers in the HD. 

We follow A05 and include in the SNIa dispersions the contri- 
bution of a peculiar velocity of 300 km/s (this dispersion depends 
on z and on cosmology) and an intrinsic dispersion of SN absolute 
magnitudes of 0.132. The latter is adjusted to obtain a SN reduced 
X?=l and coincides with the value of 0.13 ± 0.02 given in A05 for 
their concordance cosmology. 



2.2 GRBs and the L iao -E pk -T 0A5 relation 

The sample of 19 GRBs with known redshifts and with the obser- 
vational information required to establish the Lj ao = C E^Tq^s 
correlation was presented in Firmani et al. ( 2006a). The rest frame 
quantities entering in this correlation are the bolometric corrected 
(in the range of 1 — 10 4 keV in the rest frame) isotropic-equivalent 
luminosity, Liso, the spectrum peak energy, E pk , and the time 
spanned by the brightest 45% of the total lightcurve counts above 
the background, To. 45. The assumed energy range for calculating 
To. 45 is 50 — 300 keV i n the rest frame. We use the recipe pro- 
posed by iReichart et alJ 1200 ll) to pass from the observed energy 
range to the assumed rest one, and the ligth -curve time binning of 
HETE-II, 164-ms (see lFirmani et all2006j) . 

Besides the redshift, the observables required to estimate Li so , 
E pk and To. 45 are: the peak flux P, the fluence F, the spec- 
trum peak energy _E p k, and its overall shape (as described by the 
iBand et aljl993l model for most cases), and the time scale of the 
brightest 45% of the total light-curve counts above the background, 
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T0.4! • Errors on the observable parameters are appropriately propa- 
gated to the composite quantities under the assumption of no corre- 
lation among the measured errors. For the c oncordance cos mology 
the correlation gives a reduced %r = 0.7 I Firmani et al. 2006a), 
which could be signaling to an overestimate in the data uncertain- 
ties. Whilst the observable uncertainties do not improve, we can not 
attempt to estimate any intrinsic dispersion for our GRB sample. 



3 THE METHOD 

The 'circularity problem' mentioned in the Introduction in princi- 
ple should not be a problem for SNIa because the brighter-slower 
and brighter-bluer correlations could be calibrated with a low— z 
sample. However, the distance estimates to high-z SNIa improve 
if the parameters of these correlations are empirically determined 
along with the \ 2 minimization, from which the cosmological pa- 
rameters are also constrained (A05; note that multi-band photo- 
metric data are necessary to apply this technique). The situation is 
therefore, at least mathematically, similar to the circularity problem 
of GRBs. This suggests us to use our Bayesian method for improv- 
ing the SNLS cosmological constraints estimated by A05. 

The basic idea of our approach is to find the best-fitted cor- 
relation on each point Q, of the explored cosmological parameter 
space [for instance fi = (f2 m , Qa)] and estimate with such corre- 
lation the scatter x 2 (^> on the HD for any given cosmology Q. 
The conditional probability P(fi|fi) inferred from the x 2 (^>^) 
statistics provides the probability for each Q given a possible re- 
defined correlation. By defining with P'(f2) an arbitrary probabil- 
ity for each f2-defined correlation, the total probability of each f2, 
using the Bayes formalism, is given by 
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P(fi) = / P(tt\tt)P'(n)dn, 



(1) 



where the integral is extended on the available Q space. Because 
the observations give a correlation for each cosmology, P'(n) is 
actually the probability of the cosmology. Consequently such prob- 
ability is obtained putting P'(fi) = P(fi) and solving the integral 
Eq.Q. This method based on the use of an iterative Monte Carlo 
technique has been presented in Firmani et al. (2005,2006b), and 
we refer the reader to such papers for more details. 

In what follows, we apply the Bayesian method to (i) the 
SNLS SNIa sample alone, and (ii) to the combination of both the 
SNIa and GRB samples. 



4 RESULTS 

Figure Q shows the HD assuming the concordance cosmology 
(f2 m =0.28, ft A =0.72 and h=0.71) for the 117 SNIa (0.015 < z < 
1.010) reported in A05 (red symbols), as well as for the 19 GRBs 
(0.17 < z < 4.5) from our sample (blue symbols). From this 
plot one sees that GRBs are a natural extension of SNIa to high 
redshifts. The observational uncertainties for GRBs are still signif- 
icantly larger than for SNIa. The residuals of both samples with 
respect to the assumed cosmology (solid curve) are shown in the 
bottom panel of Fig.Q The average of the absolute values of the 
residuals and its uncertainty for the SNIa and GRB samples are 
0.15 ± 0.01 and 0.26 ± 0.05, respectively. 

In Fig.|2|we show the la confidence levels (CL's) for only the 
SNIa sample (green line) and for the combined SNIa+GRB sample 
(red shaded region) in the f2 m -f2A plane. The use of the Bayesian 



I I I I 



I I I I 



I I I I 



I I I I 



I I I I 




Figure 1. Top panel: The SNIa (red symbols) and GRB (blue symbols) 
Hubble diagram (HD) for a concordance cosmology. In the bottom panel 
we show the residuals of the data points minus the concordance model. 



method for analyzing the first-year SNLS SNIa dataset improves 
somewhat the cosmological constraints, especially for the large Qa 
part of the CL, as can be appreciated by comparing Fig.|2|with Fig. 
5 in A05. From Fig.|2|we also see that the inclusion of GRBs greatly 
improves the constraints given by SNIa alone. Because GRBs span 
a wide range of z's, the degeneracy between Qa and Q m is less se- 
vere fo r them than for the SNIa (see the discussion in lFirmani et alJ 
2006b). This achievement is obtained despite the small number of 
GRBs and their relatively large observational uncertainties. 

Both the SNLS SNIa sample and the GRBs sample, show 
each one that the best-fit values of S7 m and Qa are close 
to the flat-geometry case: the concordance model is actu- 
ally well inside the corr esponding la CL c onstraints (A05, 
iNesserislsT Perivolaropoulos 200ol lFirmani et alj2006bl) . Now, the 
combined set makes the constrain even more restrictive. If one 
forces f2tot= L our statistical analysis gives Q m = 0.2731q 024 ■ 
This range intersects the range of Q lrl values allowed by dynamical 
determinations (e.g., lHawkins et ai]|2003t ISchuecker etafll2003l) . 
Thus the constraints to the A cosmology parameters obtained here 
are consistent with several other independent cosmological mea- 
surements. 



4.1 Flat cosmology with alternative DE 

We now relax the assumption w = — 1, which was implicit up to 
now, and explore the possibility of w — wo, where wo is a free pa- 
rameter. The limited number of objects in our two samples and the 
current accuracies do not allow to have more than two free param- 
eters. Therefore, we fix f2 m +£7A=l, and find the CL contours in the 
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Figure 2. Constraints on the plane from the SNIa HD using our 

Bayesian approach to circumvent the circularity problem (green line) and 
from the combined SNIa+GRB HD (red line and shaded region). Both lines 
are contours at 68.3% CL's. The star shows where Xr reaches its minimum, 
while the cross indicates the concordance cosmology. The line corresponds 
to the flat geometry cosmology, the upper curve is the loitering limit be- 
tween Big Bang and no Bing Bang models 
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Figure 3. Constraints on the (O m , wo) plane for a flat cosmology with 
static DE, using the same convention of Fig.l2l 



Figure 4. Constraints on the (wo,wi) plane for a flat cosmology with dy- 
namic DE, zt=l, and f2 m =0.28, using the same convention of Fig.|2] The 
encapsulated panel shows the lcr CL for to(z)corresponding to the SNIa 
dataset alone (green plus red region) and the SNIa+GRB (red region). 



wo — plane. Figure[3]shows the la CL's in the Wq — fl m plane 
using only the SNLS SNIa sample (green line), and the combined 
SNIa+GRB sample (red line and shaded region). Again, the SNIa 
constraints obtained with our iterative Bayesian method are tighter 
than the ones obtained in A05 (compare Fig. [3] with their Fig. 5). 
When using the combined SNIa+GRB sample, we obtain a tight 
constraint on wo for reliable values of fi m . For values of f2 m in 
the range 0.236-0.286, the A model (w = —1) is consistent at the 
lcr. Assuming the prior Q, m = 0.28, we obtain u)q = — 1.055to 029' 



while for f2 m = 0.26, we obtain too = — 1.000_ 



By combin- 



ing independent cosmological probes that are sensitive to f2 m (e.g., 
lAllen et all2(X)4lEisenstein et alj2005l) with our joint GRB+SNIa 
probe, better constraints on wo could be obtained. 

In order to explore the constraints on a possible evolution o f 
w, based on the same argu ments given inlFirmani et alj J2006b), 
we use the parametrization iRaoetti et al. 2005) 



w(z) = wo + Wl 



z t + z 



(2) 



Constraints on (u>o,Wi) for flat geometry dynamical DE models 
with the priors f2 m =0.28 and Z t=l are plotted in Fig.|4| As dis- 
cussed in Firma nTet alJ l2006bl) . current observational data do not 
allow yet to determine well zt (zt = 1 cor responds to a popular 
parametrization introduced by Chev allier & Polarski 11200 H) . The 
green line and the red shaded region represent the la CL's for the 
SNIa dataset alone and for the combined SNIa+GRB sample, re- 
spectively. The encapsulated panel shows the corresponding la CL 
for the evolution of w(z) using the SNIa dataset alone (green plus 
red region) and the SNIa+GRB datasets (red shaded region). Again 
the A case (wo=— 1 and wi=0, which reduces to the concordance 
model because of the assumption that fl m = 0.28), is within the la 
CL's for both SNIa and SNIa+GRB constraints. Our results allow 
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at the la CL for models that avoid crossing the phantom dividing 
line (w [z] = —1). Notice that although Eq. describes the evo- 
lution of w up to any arbitrary large z once its parameters were 
determined, the changes on w with z suggested by the observa- 
tional constraints are formally valid only within the redshift range 
of the observational data, z < 4.5 in our case. 



5 CONCLUSIONS 

We have combined a sample of potential standard candles that in- 
cludes 1 15 SNIa of the SNLS dataset (A05) and 19 long GRBs. The 
latter were "standarized" on the basis of a ti ght correlation among 
prompt 7— ray properties alone iFirmani et alj|20 06a). Exploiting 
some similarities between the energetic callibrations of the SNLS 
SNe and the GRBs, we use for both populations of objects the same 
method to constrain cosmological parameters, namely a Bayesian 
approach described in Firmani et al. (2005,2006b). GRBs may be 
conceived as the natural extension of SNIa to large z's in the HD. 
The advantage of this extension by using GRBs is that the heavy 
degeneracy in the space of the cosmological parameters due to a 
narrow and low z range (as is the case of SNIa) is remarkably re- 
duced. Our main results are as follow: 

• The cosmological constraints obtained with the Bayesian 
method for the SNLS sample alone are in agreement with those 
reported in A05. However, with the Bayesian method we obtain 
somewhat tighter constraints on fi m , JIa and wo than with the \ 2 
minimization method used in A05. The values of Sl m and Ha of 
the best fit (using SNIa only) are in good agreement with the flat- 
geometry case. Moreover, the best fitting values of fi m and Oa 
do not disa gree with those obta ined by using other cosmological 
probes (e.g., Spereel et al. 2006), re-affirming the reliability of the 
cosmological concordance model. 

• We have presented in a previous paper iFirmani et all20 06b) 
cosmological constraints derived by using GRBs only, whose ener- 
getics is standardized with the Ljso--Epk-To.45 relation. We have 
shown here that the combined SNIa+GRB sample is able to largely 
reduce the allowed region of the parameter space with respect to the 
case where a single population is used. The resulting values of the 
best fits, once again, are in agreement with the concordance ACDM 
model. 

• As a consistency check, we have explored DE models with 
w = u>o =const (i.e. relaxing the assumption of w = —1), but 
assuming flat geometry. Furthermore, for completeness, we have 
allowed also w to change with z according to a parametric law and 
assuming fi m =0.28. In both cases we find that w = — 1 =const. is 
within the la CL from the best fits. 

Finally, we emphasize the next two general conclusions: 

1) GRBs and SNIa as standard candles should not be con- 
sidered as competing cosmological probes but as complementary 
methods. Besides, both GRBs and high-z SNIa suffer from the 
circularity problem concerning their calibration in such a way the 
same Bayesian method can be applied for both samples. 

2) According to our results there is no need for DE "exotic" 
equations of state. The flat Friedmann-Robertson-Walker-Lemaitre 
A cosmology is fully consistent with the HD constructed for the 
joint sample of SNIa and GRBs up to z = 4.5. Similar conclu- 
sions were obtained recently on the basis of other high-precision 
cosmological probes 1 Sper gel et all2006l) . 
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